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Abstract 

Spectral stability of multi-hump vector solitons in the Hamiltonian system of coupled nonlinear 
Schrodinger (NLS) equations is investigated both analytically and numerically. Using the closure theo- 
rem for the negative index of the linearized Hamiltonian, we classify all possible bifurcations of unstable 
eigenvalues in the systems of coupled NLS equations with cubic and saturable nonlinearities. We also 
determine the eigenvalue spectrum numerically by the shooting method. In case of cubic nonlinearities, 
all multi-hump vector solitons in the non-integrable model are found to be linearly unstable. In case of 
saturable nonlinearities, stable multi-hump vector solitons are found in certain parameter regions, and 
some errors in the literature are corrected. 

1 Introduction 

The coupled NLS equations have wide applications in the modeling of physical processes. For instance, such 
equations with the cubic nonlinearity govern the nonlinear interaction of two wave packets 0] and optical 
pulse propagation in birefringent fibers or wavelength-division-multiplexed optical systems 115) . 
Similar equations with the saturable nonlinearity describe the propagation of several mutually-incoherent 
laser beams in biased photorefractive crystals [111 HI)] . Various types of vector solitons including single- 
hump and multi-hump ones have been known to exist in these coupled NLS equations [21 El El HOI 1111 1141 
I16L I25L 1^ 1341 and they have been observed in photorefractive crystals as well |^ 171 124| . 

Linear stability of vector solitons in the coupled NLS equations is an important issue. Fundamental single- 
hump vector solitons are known to be stable |2(il 1291 r?7] . Stability of multi-hump vector solitons (which 
have one or more nodal points in one or more components) is more subtle. For the cubic nonlinearity, it 
was conjectured in |36| based on the numerical evidence that multi-hump vector solitons were all linearly 
unstable. If the multi-hump solitons are pieced together by a few fundamental vector solitons, then their 
linear instability has been proven both analytically and numerically in j^Hl QH] ■ The linear instability for 
other types of multi-hump vector solitons has not been proven yet. For the saturable nonlinearity, multi- 
hump solitons have been shown to be stable in certain parameter regions j^SlEHli but the origins of their 
stability and instability have not yet been fully analyzed. 

^From a broader point of view, the theory of linear stability of vector solitons in coupled NLS equations was 
recently developed with the use of the closure theorem for the negative index of the linearized Hamiltonian 
UHl EH- However, there are not many applications of the general theory to particular bifurcations of 
unstable eigenvalues jEll^i because the general theory excludes non-generic bifurcations. It is desirable 
to further develop a perturbation theory to the eigenvalue bifurcations, so that the origin of instability 
becomes more apparent in the context of the closure theorem. 



In this paper, we investigate the hnear stabihty of multi-hump vector soHtons in the general Hamiltonian 
system of coupled NLS equations both analytically and numerically. Using the closure theorem for the 
negative index of the linearized Hamiltonian as well as the perturbation technique, we classify all possible 
bifurcations of unstable eigenvalues in two physical models with cubic or saturable nonlinearities. In 
the first model, we show that multi-hump vector solitons near points of local bifurcations are always 
linearly unstable, in agreement with numerical results in |36| . In the second model, the situation is more 
complicated. Our results show that the 1st family of multi-hump vector solitons is indeed linearly stable near 
the local-bifurcation boundary, in agreement with numerical results in |26j . However, for the 2nd family, we 
discovered a new oscillatory instability near the local-bifurcation boundary, which was missed in 26 . Due 
to this oscillatory instability, the stability region of vector solitons for the 2nd family is drastically reduced 
from that reported in . Numerically, we track the unstable eigenvalues of multi-hump solitons and reveal 
various scenarios of eigenvalue bifurcations away from the local-bifurcation boundaries. We also map out 
the correct stability regions of multi-hump vector solitons in the entire parameter space. Furthermore, the 
number of numerically-obtained unstable eigenvalues agrees completely with that predicted by the negative 
index of the linearized Hamiltonian. 

Our paper is structured as follows. The main formalism and the closure theorem for the negative index of 
the linearized Hamiltonian are described in Section 2. Analysis of unstable eigenvalues in the coupled NLS 
equations with cubic and saturable nonlinearities is developed in Sections 3 and 4, respectively. Section 5 
summarizes our results and open questions. Appendix A reviews bifurcations of unstable eigenvalues by 
the perturbation method. 



2 Main formalism 

We consider a general Hamiltonian system of coupled NLS equations in the form: 

where z S K+, a; G M, V« G C, dn G M, and U = C/(|-0iP, IV'JVP) G K. We assume that U{0) = U'{0) = 0, 
and dn > for all n. In optical fibers (photorefractive crystals), the function ^„(z,a;) is the envelope 
amplitude of the n*-^ channel (beam), z is the propagation distance along the fiber (waveguide), and x is 
the retarded time (the transverse coordinate) ^2 El 1161 EH] ■ 

Following the recent work in [22], we study the linear stability of vector solitons: 

V'„(2,a;) = $„(x)e*''"^ (2.2) 

where $„ : R ^ M, and /3„ > for all n. We assume that none of the components $„(a;) vanish identically 
on cc G M. Linearization of the coupled NLS equations 1)2. l|l follows from the expansion: 

ipniz,x) = |$„(a;) + [un{x) +iwn{x)]e^'' + [un{x) + iwn{x)] e^^|e''^"^ (2.3) 

where ||u„|j , ||?x;„|j <^ 1, and the overline denotes the complex conjugation. The linearized equations for 
(u„, Wn) are the following non-self-adjoint problem in L^(]R, C^^): 

£iu = —Aw, Co'w = Au, (2.4) 

where A G C is an eigenvalue, (u, w)-^ : M C^^ is the eigenvector, and Co and £i are the matrix 
Schrodinger operators with elements: 

f dU \ 
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f dp dU \ d'^U 

(Ak™ = (^-^"^+/?"-5^j^"--2^^^ci>„$„. (2.6) 

Since the eigenvalue problem H2.4|l is a linearization of the Hamiltonian system, the values of A occur as 
pairs of real or purely imaginary eigenvalues, or as quadruplets of complex eigenvalues. Eigenvalues with 
Re(A) > lead to spectral instability of vector solitons H2.2|l . We denote the number of eigenvalues in the 
first open quadrant as A^comp, the number of positive real eigenvalues as iVrcai, and the number of purely 
imaginary eigenvalues with positive Im(A) as iVimag- The continuous spectrum has N branches, located at 
the positive imaginary axis for Im(A) > Pn, n — I, N . Zero eigenvalue A = has geometric multiplicity 
of at least (iV + 1) and algebraic multiplicity of at least {2N + 2), in the assumption that none of the 
components $„(a;) vanishes identically on x € M \'27\ . 

Furthermore, we denote the number of negative and zero eigenvalues of operators £o,i in £^(M,C^) as 
n(£o,i) and z{Cq.i), respectively. We also assume that the solution $ri(x) depends smoothly on ...,f3N) 
in an open non-empty set of and introduce the Hessian matrix U with elements: 

l^n,m = ^r, 1 (2-7) 

where Qn — QniPi, ■■•,/3jv) = /r '^'n'^^;. We denote the number of positive and zero eigenvalues of matrix 
U as p{U) and z{U), respectively. Finally, we introduce the linearized Hamiltonian ("energy") of the 
eigenvalues A in H^{R, C^^): 

/i[u,w] = (u,Au) + (w,/:ow), (2.8) 

where (•, •) is the standard inner product in L^(]R, C^^). The negative index of the linearized Hamiltonian 
is the number of negative eigenvalues of Ci and Cq in L^(R, C^). 

Several assumptions are imposed on the linearized problem (|2.4|l in a general case |27| : 

(i) z(£i) = 1, ziCo) = N; 

(ii) z{U) - 0; 

(iii) no eigenvalues A e iM exist with h[u,w] = 0; 

(iv) no embedded eigenvalues A e exist with |Im(A)| > /Jmin, where Pmhi = niin(/3i, I3n)- 

Closure Theorem |27| Assume that (i)-(iv) be satisfied. Let N-^^^ be the number of eigenvalues A G iR-|- 
with /i[u, w] < 0. Then 

(i) Aureal + 2A^comp + 2A^;,,g = n(£i)-p(Zi)+n(£o), (2.9) 

(ii) Nrc.i>\niCi)-p{U)-n{Co)\, (2.10) 
(in) N,,^p <mm (n{Cn),n{C^)~p{U)), (2.11) 

such that 

\n{C,)-p{U) - n(£o)| < A^unst < n(£i) - p(Z^) + n(£o), (2.12) 

where iVunst = Apical + 2A'comp is the total number of unstable eigenvalues in the problem X'li.J^ . 

This theorem was originally proved for the coupled NLS equations H2.1|l in one dimension (27| and then 
generalized to a three-dimensional NLS equation [5] and to an abstract Hamiltonian dynamical system 
|18j . It allows us to analytically trace unstable eigenvalues under parameter continuations, starting with 
the particular limits, where all eigenvalues A of negative energy ft,[u, w] are known. Examples of such 
parameter continuation are recently reported in [171 141| in the context of the coupled NLS equations. 

Bifurcations of unstable eigenvalues may occur in the linearized problem H2.4|l . when operators £i and £o 
change according to a continuous deformation and one of the assumptions (i)-(iv) of the Closure Theorem 
fails. Bifurcations are reviewed in Appendix A. In what follows, we apply parameter continuation and 
bifurcation analysis to the system of coupled NLS equations (|2.1|) with cubic and saturable nonlinearities. 
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3 The coupled cubic NLS equations 



We consider the system of coupled cubic NLS equations [H 1151 1^ : 

#2z+l/'2.. + (xlV'l|' + IV'2|')V'2 = 0, (3.1) 

where x > 0. The system is a particular example of H2.1|l with N — 2, di — d2 — 1, and 

U=lHi\' + x\i^i\^\H' + l\H^- (3.2) 

The system has a countable infinite set of families of vector solitons 4>(a;) = ($i, $2)"^, classified by 
different nodal index i = {ii,i2)'^, where i„ is the number of zeros of <i>„(x) on .x G M ^1 114L 15^ . We 
consider here families of vector solitons with nodal index i = (0, n)"^, n G N, which are locally close to 
the NLS soliton, *NLs(a:) = ($(°),0)^, where $(o)(a;) = sech(VAa;). We let /3i = 1 and /Sa = /? for 

convenience and introduce scalar Schrodinger operators: 

Lo = -^ + l-2sech^(x), (3.3) 

Li = __ + l_6sech2(x), (3.4) 

Ls = -^ + A2(x)-2xsech2(a;), (3.5) 

where 

Mx) = ^ ■ (3-6) 

The scalar operators Lo,Li,Ls define the matrix operators Cq and £1 at e = 0: 

Co = diag(Lo, L^), Ci = diag(ii, L^). 
We define the perturbation series expansions of vector solitons $ = ($1, $2)^: 

$i(x) = $(°)(x)+e2$(2)(a;) + 0(e*), 

$2(2;) = e<^>^'^\x)+e^<^^^\x) + 0{e^), (3.7) 

and 

(3^Xlix) + e'Cn{x) + 0{e'). (3.8) 
Corrections of the perturbation series (|3.7|) - (|3.8|) satisfy the linear equations: 

= 0, (3.9) 
^ x*'"^ ($*'^)', (3.10) 
L,$(3) = -C„(x)$(') +2x$(°)$W$('^ + ($(^^)'. (3.11) 

The problem H3.9|l has a decaying solution <i>*^^^ = ^n\x) (see |3B]). When n = and x > 0, the 
solution — sech'*(a;), s = Ao(x) is a ground state. When n > and X > Xn = "-(" + l)/2, the solution 
^nHx) is an excited state with exactly n nodes on a; £ K. The problem (|3.1U|I also has a decaying solution 
$'^^(x), since the right-hand-side orthogonal to the kernel of the operator Ci, which is 

$'"^'(2;). By the Fredholm Alternative Theorem, the problem (|3.11f) has a decaying solution if and only if 
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the right-hand-side is orthogonal to the kernel of Cg, which is (^n\x). The orthogonality condition defines 
the parameter C„(x) in the form: 



)^2 



Cnix) = \ n) (1). —■ (3-12) 

The condition C„ (x) 7^ gives the sufficient condition of continuation of the perturbation series expansions 
(|3.7|I - H3.8|I . Thus, for x > x„ and C„(x) ^ 0, there exists some i?„ > 0, such that the n-th family of 
vector solitons '^{x) = ($1, $2)^ with the nodal index i — (0, n)'^ bifurcates from $nls — 0)^ in the 

one-sided domain Bn'- 

6„ = {/3:0< |/3-A2(x)| <i?„, sign (/3 - (x)) = sign(C„(x))} ■ (3.13) 

These results for the first three families n = 0, 1,2 were analytically obtained and numerically verified in 
We investigate stability of the n-th family of vector solitons in the one-sided domain Bn below. 



3.1 Analytical results 

We trace unstable eigenvalues using the Closure Theorem. We consider a generic case C„(x) 7^ in the 
one-sided open domain Bn and show that the left-hand and right-hand sides of the closure relation H2.9|l 
are equal to 2n for small e > 0. 

Operator Lq in H3.3|) has one bound state for zero eigenvalue, operator Li in H3.4|) has two bound states for 
negative and zero eigenvalues, and operator in H3.5|l has in + 1) bound states with n negative and one 
zero eigenvalues. Therefore, at e = 0, we have n(£o) = + n = n, z{Ca) = 1 + 1 = 2, n{Ci) = 1 + n and 
z{Ci) = 1 + 1 = 2. It follows from Sturm Nodal Theorem that 



n{Co) = n, z{Co) = 2, 



Ve > 0. 



Since z{£i) = 2 > 1, we have the bifurcation case z{Ci) > 1 for e = (see Appendix A.l). It is however a 
degenerate bifurcation case, since it occurs on the boundary of the existence domain Bn, such that f3 S dBn- 
We trace the zero eigenvalue of £1 for e 7^ by the regular perturbation series. 



u(x) 








' M(i)(a;) ' 









+ e 










and 



Corrections of the perturbation series (|3.14(l satisfy a set of linear non-homogeneous equations: 

Liii(i) = 2x'i>(°)('i>i^))', 

L,u(^^ = (A2-a.(xM^)+2x<i>(°)<i>«(z.«+<i>(^))+ 3 



(3.14) 

(3.15) 

(3.16) 
(3.17) 



It follows from (|3.1Q(I and (|3.16(l that m^^^ = 2<i>(^). By the Frcdholm Alternative Theorem, decaying 
solutions of H3.17|) exist if and only if the right-hand-side of (|3.17|) is orthogonal to ^n\x). Using (|3.12|) . 
we find that A2 = — 2C„(x). Therefore, we have: 



n(£i) = l + e(C„(x)) + n, z(£i) = 1, 



e > 0, 
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where 0(z) is the Heaviside step- function. We trace the zero eigenvakie of U from (|2.7(l and (|3.7|l : 



= 

U2,l = 
1^2,2 = 

It follows from (|^^ that 

and, due to the symmetry of U, 



dQi 



d(3i 
dQi 



= 2 + 2($(°\$(2)' 



+ 0(6^), 



dP2 
dQ2 



2) 



dP2 



/3i = l 



/3l = l 



/3i = l 



^W2 

df3i 



($(1) $(!))_ 



2($(o),$( 



($(1) $(1) 



/3i = l 



o(.^), 



/3i = l 



+ 0(6^). 



Pl = l 



5/32 
det(Z^) 



/3i = l 



Cnix) 



C„(x) 



(3.18) 



(3.19) 



Therefore, we have: 



e > 0. 



We conclude that the bifurcation case z{Ci) > 1 on the boundary of the existence domain (3 G 9i3„ does not 
result in bifurcation of any eigenvalue A of the stability problem H2.4|l . such that n{Ci) —p{U) -|-ri,(£o) = 2n 
is valid everywhere in /? g i3„ U dBn- It follows from the Closure Theorem that the ground state with n — Q 
is spectrally stable in /? e H„, while the n-tli excited state with n > \ may have at most A^unst unstable 
eigenvalues, where < Nu_nBt < 2?t,. We show that A^unst — 2A^comp = 2n in /3 G Bn in a generic case. 

At e = 0, the stability problem H2.4|l can be decoupled as follows: 



-\w^ 



LqWi — Ami 



and 



Ls(u2 ± iW2) = ±iX{u2 ± iW2) 



(3.20) 



(3.21) 



The first problem H3.20|l has the continuous spectrum for Re(A) = and |Ini(A)| > 1 and the zero eigenvalue 
A = of algebraic multiplicity 4 and geometric multiplicity 2. The second problem H3.21|l has the continuous 
spectrum for Re(A) — and |Ini(A)| > A^(x), zero eigenvalue A = of geometric and algebraic multiphcity 
2, and 2n isolated eigenvalues in the points A = ±i (A^ — A^), where fc = 0, 1, n—1. It follows from l|3.6|l 
that for X > Xn- 



4 - = {n~k) [2A„(x) + {n- k)] > {n - kf > 1, 



< fc < 



(3.22) 



Therefore, 2n isolated eigenvalues of the problem H3.21|l are embedded in the continuous spectrum of the 
problem (13.20(1 . These embedded eigenvalues have negative energy /i[u,w], since at e = 0: 



0<k < 



(3.23) 



where — (0,$^ ) and Wfc = (0, =F**J?fc ) at e = 0. By Appendix A. 4, all 2n embedded eigenvalues 
of negative energy /i[u, w] bifurcate in a general case of non-zero T, see Eq. I|A.29|) . to complex unstable 



eigenvalues A G C, Re(A) > for e 7^ 0, such that A^ioai + 2Afcomp + 2N.] 



2N,. 



2n in /? G Bn 
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3.2 Numerical results 



For e ^ 0, the linearized problem (|2.4|) satisfies the assumptions of the Closure Theorem. Therefore, 
all unstable eigenvalues A'^Linst = SA'comp — 2n are structurally stable for larger values of e, until new 
bifurcations occur in the parameter continuations. We study numerically locations of unstable eigenvalues 
in the linearized problem 12.4|l related to the vector solitons $ = {^1,^2)^ with nodal index i = (0,n)^, 
n ~ 1,2. Our numerical algorithm is based on the shooting technique in the complex A-plane. We also 
determine the indices n{CQ), n{Ci) and p{U) by a numerics-assisted procedure as described in |39|. and 
relate them to the number of unstable eigenvalues by using the closure relation H2.9|l . 

Figure 1 shows the Ist-family of multi-hump vector solitons with the correspondence: u = v = $2(x), 

and LU = y/j3. For uj < 1, this family exists between Xi(/3) < x < X2{l3), where x — X2(/3) is the local 
bifurcation boundary, and X = Xi iP) is the nonlocal bifurcation boundary jS] . Hence the one-sided domain 
P G Bi is located to the left of the local bifurcation curve, and sign(d) = 1 in When parameter 

UJ = 0.6 is fixed, we readily find that xi ~ 0.28 and X2 = 2.08. When x moves from X2 to xi, the distance 
between the two pulses in the v component grows. It diverges to infinity at the nonlocal bifurcation 
boundary X2 (near point a in Fig. 1). 

Figure 2 shows unstable eigenvalues of the linearized problem (|2.4f) for uj = 0.6 and Xi < X < X2- In the 
domain f3 G Bi, there is a pair of unstable complex eigenvalues A ~ Re(f72) ± iIm((T2), which bifurcate from 
the embedded eigenvalues A = ±i(AQ — A^), see Eq. (|3.22l) . 

When X — 1^, the complex eigenvalues A = Re((T2) ± iIm(CT2) approach the imaginary axis and become 
embedded eigenvalues A — ±i{l — uj'^). The case x — ^ corresponds to the integrable Manakov system, 
when the linearized problem (|2.4I) has the following exact solution: 



uo= ^ ^o = T^{ J ], \ = ±i{l~to^). (3.24) 





This exact solution is generated by the additional polarizational-rotation symmetry in the potential function 
U — [/(IV'iP + |'02p) at X = 1- Embedded eigenvalues A = ±i(l — ui'^) have negative energy h[u,w], since 

(uo,/:iUo) - (wo,/:oWo) =-(l-cc;2)(($i,$i)- ($2,$2)) <0, (3.25) 
where the last inequality is confirmed numerically from integration of the exact solutions |34) : 



$ ( ) \/l — coshcja; 

cosh X cosh LUX — uj sinh x sinh wx ' 



$ ( ) wyT— c^sinhx ^j-j 

cosh x cosh ujx — uj sinh x sinh ujx 

in the entire domain of existence: < w < 1. Since embedded eigenvalues A — ±i(l — uj"^) at x = 1 have 
negative energy /i[u,w], they bifurcate to the complex plane for x 7^ 1 when the polarizational symmetry 
is destroyed. This is indeed the case as shown in Fig. 2, in agreement with Appendix A. 4. 

There is an additional instability bifurcation at x = 1- This bifurcation comes about because at this x 
value, the Ist-family of vector solitons # — (<&i,$2)"^ can be generalized to asymmetric solutions with an 
additional free parameter [51 114L OS] . As a result, the derivative of the asymmetric vector solitons with 
respect to the free parameter is an eigenvector in the kernel of the operator Ci, such that z{Ci) = 2 at 
X = 1. When x 7^ 1> the integrability of the Manakov system is destroyed, and a pair of real or purely 
imaginary eigenvalues is generated, in agreement with Appendix A.l. Indeed, Fig. 2 shows a pair of purely 
imaginary eigenvalues A = ztiai for x > Ij which merge to the end points A = ztiuj'^ of the continuous 
spectrum at x = 1.185, and a pair of real eigenvalues A = ±ai for x < 1- 
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Now we relate results of Fig. 2 to the closure relation H2.9|l . For this purpose, we have determined the 
indices n{£i) and p{U) by the numerics-assisted procedure in [351 (we note that n(£o) = 1 everywhere in 
the existence domain of the Ist-family of vector solitons). For to = 0.6, we have found numerically that 

n{C,) = I ;J' ^1 < ^ < 1' p{u) ^ 2, for all Xi < X < X2, (3.28) 
3, l<x<X2, 



such that 



n{C^) + niCo)-p(U)^\ I' < ^ < 1' (3.29) 

2, l<x<X2- 



On the other hand. Fig. 2 shows that 



iVcomp = 1, for all Xl<X< X2, iVrcal = <! I' ^_ ^ ^_ ^ (3.30) 



0, 1< X < X2 



and the closure relation H2.9|) is thus satisfied. 



Figures 3 and 4 show similar results for the 2nd family of multi-hump vector solitons. When parameter 
oj — 0.6 is fixed, the local bifurcation boundary is X2 — 4.68 and the nonlocal bifurcation boundary is 
Xi = 1.68. Again, the one-sided domain f3 E B2 is located to the left of the local bifurcation boundary, such 
that sign(C2) = 1. However, such solitons exist on both sides of the nonlocal bifurcation boundary x = Xi- 
As X nioves leftward from X = X21 it first crosses the nonlocal bifurcation boundary x — Xi? then turns 
around, and approaches the nonlocal bifurcation boundary from the left side (see point a in Fig. 3). This 
behavior has been explained both analytically and numerically in 

Using the same shooting algorithm, we have obtained the unstable eigenvalues of the linearized problem 
l|2.4|l and displayed them in Fig. 4. In the domain f3 E B2, there exist two pairs of unstable complex 
eigenvalues, such that the pair A = Re(cr4) ± iIm((T4) bifurcates from the pair of embedded eigenvalues 
A = ±«(Ao — A2), while the pair A — Re((T3) ± ilm(cr3) bifurcates from the pair of embedded eigenvalues 
A = ±«(A^ — A|). Fig. 4 also shows that the pair A = Re(tT3) ± zlm(cr3) approach the imaginary axis and 
become a pair of embedded eigenvalues at x = Xa ~ 2.49, but then reappear as a pair of complex unstable 
eigenvalues for x < Xa, in agreement with Appendix A. 4. 

There are two more instability bifurcations in Fig. 4. At x = Xfc ~ 2.44, the zero eigenvalue bifurcates into 
a pair of imaginary eigenvalues A = ±i(T2 for x > Xfci which then merge into the end points A = ztiuj'^ of the 
continuous spectrum at x ~ 2.72. When x < Xb, this zero eigenvalue bifurcates into a pair of real unstable 
eigenvalues A = ±172. At yet another point X — Xc — 2.17, the zero eigenvalue bifurcates into a pair of 
imaginary eigenvalues A — zLiai for x > Xci which merge into the end points A — zLiuj'^ of the continuous 
spectrum at x ~ 3.01. When x < Xc, this zero eigenvalue bifurcates into a pair of real unstable eigenvalues 

A = ±<Tl. 

To relate numerical results of Fig. 4 to the closure relation H2.9|l . we have again determined the indices 
n{Ci) and p{U) by the numerical algorithm, while u^Cq) ~ 2 everywhere in the existence domain of the 
2nd family of vector solitons. For lu — 0.6, we have found numerically that 

niC) = I ^ < piU) = I ^' ^ <>^- (3.31) 

[ 4, Xh < X < X2, 2, Xc < X < X2, 



such that 



6, Xi < X < Xc, 

n{Ci)+niCo)~p{U) = { 5, Xc < X < Xb, (3.32) 

4, Xb < X < X2- 
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On the other hand, Fig. 4 shows that 

{2, Xi < X < Xc, 
1, Xc < X < Xb, (3-33) 
0, Xfc < X < X2- 

Hence the closure relation (|2.9|l is satisfied. In particular, the instability bifurcation a.t x — Xb is due to a 
jump in the index n(£i) (similar to the 1st family), while the instability bifurcation at x = Xc is due to a 
jump in the index p{U). These two bifurcations occur in agreement with Appendices A.l and A. 2. 

Although our results were obtained here for a particular value lo = 0.6, we expect that similar results hold 
for other values of w, when < < 1. We conclude that the 1st and 2nd families of multi-hump vector 
solitons in the coupled cubic NLS equations are all linearly unstable (except for the integrable Manakov 
system x = 1, where the 1st family is neutrally stable). 

It is remarkable that the main features of instability bifurcations for the Ist-family of multi-hump vector 
solitons are repeated for the 2nd-family of vector solitons, irrelevant whether the coupled NLS equations are 
integrable or not. This allows us to conjecture that a similar pattern of unstable eigenvalues persists for a 
general n-th family of multi-hump vector solitons, with more unstable eigenvalues and additional instability 
bifurcations appearing as n increases. 



4 Two coupled saturable NLS equations 

We consider the system of two coupled saturable NLS equations ^] El '■ 



IV'1P + IV'2P 
l + s(|VlP+|V'2 
IV'lP + IV'2p 



i^2z+^2xx + ^ , /I / 12 I I / I2\ ^2 = 0, (4.1) 

where s > 0. This system is a particular example of H2.1|l with N ^ 2, di = d2 = 1, and 

U=\ (l^iP + 1^21') - ^log (1 -f 5(|^ip -t- |V^2n) . (4.2) 

We consider again the n-th family of vector solitons $ = ($i,<&2)"^ with nodal index i = (0, n)-'", n e N 
and convenient parametrization (3i = \ and (32 = (3- The n-th family bifurcates from the scalar solution 
* = ($0,0)^, where ^o{x) satisfies the ODE: 

*o-*o + Y^ = 0. (4.3) 

The local bifurcation occurs at P — (3n{s), when there exists a n-nodal bound state $ra(x) in the linear 
eigenvalue problem: 

<-/3A + T^^-0- (4-4) 



Vector solitons in the n-th family disappear at the nonlocal bifurcation boundary. The domain of existence 
for the first three families n = 0, 1, 2 has been obtained numerically in [SH^Hl- We trace analytically unstable 
eigenvalues and show that Apical + 2A^comp + 2A^j^,ag ~ 2n for small |/3 — /9„(s)| ^1. At /3 = Pnis), the 
stability problem 1)2. 4|l can be decoupled as follows: 

LiUi — —\wi, LqWi — Xui (4-5) 
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and 



Ls{u2 ± iW2) = ±iX{u2 ± iW2), 



(4.6) 



where 




dx^ 1 




(4.7) 




(4.8) 



(4.9) 



The first problem (|4.5(l is the linearized stability problem in the scalar saturable NLS equation (|4.1|l for 
^o{x). Based on numerical data in [221123 j we assume that the bound state ^o{x) is spectrally stable in 
the scalar saturable NLS equation and the problem H4.5|l does not have any eigenvalues of negative energy 
h[u, w]. It has the continuous spectrum at Re(A) = and |Im(A)| > 1, the zero eigenvalue A = of algebraic 
multiplicity 4 and geometric multiplicity 2, and possibly isolated eigenvalues A e iR of positive energy. 

The second problem H4.6|l has the continuous spectrum at Re(A) = and |Im(A)| > (3n{s), zero eigenvalue 
A = of geometric and algebraic multiplicity 2, and 2n isolated eigenvalues A = ±i[Pk{s) — Pn{s)], where 
k — 0,l,...,n— 1. By the Sturm Nodal Theorem, eigenvalues (3k{s) (k = 0,...,n) are ordered in the 
decreasing order and are characterized by eigenfunctions $fe(a;) with k nodes, such that the ground state 
$0(2^) corresponds to /3o{s) = 1 and the n-th excited state "I>„(a;) corresponds to f3n{s). The 2n eigenvalues 
have negative energy ft,[u, w], since 



where = (0, $fc)^ and Wfe = (0, =F*$fc)^ at ^ = /3„(s). Therefore, ^N-^^^^^ = 2n at ^ = /3„(s). Since the 
zero eigenvalue has the generic algebraic multiplicity six, no negative eigenvalues of Ci and Cq arise from 
the zero eigenvalue, such that we have A'^roai + 2iVcomp + "^N^^^ = 2n for \(3 — /J„(s)| <C 1 by continuity of 
the negative index n{h). The ground state with n — is therefore spectrally stable in the existence domain 
(which is /3 = 1, s > for n = 0). 

We show that the nodal bound state with n > 1 may have at most A^unst unstable eigenvalues, where 
< iVunst < 2n~ 2. The number of unstable eigenvalues is reduced by two, since there are two eigenvalues 
A = ±i(l — /3) of negative energy, which exist for all vector solitons with n > 1 due to the polarizational- 
rotation symmetry in the potential function U = C/(|^ip + |'02p)- This pair of eigenvalues is similar to 
the one which occurs in the coupled NLS equations with x = 1; such that the stability problem H2.4|l 
has exactly the same solution (|3.24|1 for /3 = oj"^. These eigenvalues have negative energy /i[u,w] due to 
H3.25|l . where the inequality remains true in the entire existence domain, as follows from numerical data in 
j26| . If /3 < 1/2, these eigenvalues A — ±i{l — (3) are embedded in the continuous spectrum of the problem 
H2.4|l . but never bifurcate off the continuous spectrum as the parameters vary. Therefore, 2A'j~ > 2 and 
< iVunst < 2n — 2 for n > 1. As a result, the 1st family of vector solitons is spectrally stable when (3 is 
near the local bifurcation boundary (3 = f3n(s). 

When (3 — /?2(s), the eigenvalues A ~ zki [(3i{s) ~ /32{s)] are embedded into the continuous spectrum if 
(3i{s) > 2/32(3), and are isolated from the continuous spectrum if /3i(s) < 2/32 (s). In the first case, 
the embedded eigenvalues bifurcate generally to complex unstable eigenvalues for < |/3 — /32(s)| <C 1, 
according to Appendix A. 4. In the second case, the isolated eigenvalues do not bifurcate to complex 
unstable eigenvalues for < \P — /32{s)\ <C 1. Since these eigenvalues have negative energy [see Eq. H4.1()|l ]. 
the number of unstable eigenvalues A^unst is then zero. In other words, vector solitons of the 2nd family 
near the local bifurcation boundary with /3i(s) < 2/32 (s) are spectrally stable. 



[f3k{s) - (3n{s)] <0, 



< fc < n. 



(4.10) 
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We compare the above theoretical results with numerical results in where unstable eigenvalues in the 
linearized problem H2.4|l were obtained for the 1st and 2nd families [see Figs. 2 and 3 in where (A,/3) 
correspond to our parameters (/3, tr)]. It was shown in that the Ist-family of vector solitons is stable 
in the domain /3i(s) < (3 < P]^\-\^{s)^ in agreement with our prediction of linear stability near the first local 
bifurcation boundary (3 = j3i{s). At /? = I3[^^^{s), the bifurcation z{lA) — 1 occurs, which generates a pair 
of real eigenvalues for (3 > /3stab('') ^^'^ ^ P^^^ '^^ imaginary eigenvalues for P < Pstlhi^)^ agreement with 
Appendix A. 2. Therefore, 

n{Ci) + niCo) - piU) = | 

and 

^imag = 2, for aU /3, iV,eal 

For the 2nd family, the eigenvalues A = ±i[/3i(s) — /32(s)] are embedded when 0.646 < s < 0.857 and 
isolated when s > 0.857. According to our prediction, the embedded eigenvalues should bifurcate to the 
complex unstable eigenvalues for < |/3 — P2{s)\ <C 1. However, it was claimed in that vector solitons 
in the 2nd family were stable near the local bifurcation boundary < |/3 — /32(s)| ^ 1 for all values of s. 
In order to resolve this discrepancy, we have numerically determined the eigenvalue spectrum for vector 
solitons in the 2nd family by the shooting method. Figures 5 and 6 present numerical results for /3 = 0.16 
and P = 0.49, respectively. 

For /3 — 0.16, the local bifurcation boundary of the 2nd family occurs at s = 0.785 and the eigenvalues 
\ — ±i [/3i(s) — /?2(s)] are embedded, such that a pair of unstable complex eigenvalues A — Re(a2)±ilm(cr2) 
bifurcates for s < 0.785. These complex eigenvalues indicate the oscillatory instability of vector solitons 
near the local bifurcation boundary |/3 — (32{s)\ ^ 1 for f3 — 0.16. Thus the claim in 26 on the stability of 
2nd- family vector solitons anywhere near the local bifurcation boundary is incorrect. Note that the complex 
eigenvalues A = Re((T2) i iIm((T2) persist throughout the entire existence domain of the 2nd family, which 
is 0.513 <s < 0.785 at /3 = 0.16. 

Furthermore, a pair of purely imaginary eigenvalues A — ±i|cri| bifurcate from the end points A = ±ip of the 
continuous spectrum at s = 0.763, merge into the origin at s = 0.696, and then bifurcate into a pair of real 
eigenvalues A — ±(Ti for s < 0.696. This exponential instability induced by the real eigenvalue A = ai has 
been reported in j26|. Since it was claimed in |26| that z{iU) = in the entire existence domain of the 2nd 
family of vector solitons, we conclude that the instability bifurcation at s = 0.696 falls into the scenario of 
Appendix A.l with z{Ci) = 2. In the interval 0.513 < s < 0.696, the oscillatory instability is overshadowed 
by the exponential instability from the eigenvalue X = ai. However, in the interval 0.696 < s < 0.785, 
this oscillatory instability is the only instability experienced by vector solitons. Since Re((T2) is less than 
0.011 in the interval 0.696 < s < 0.785, it may explain why this oscillatory instability was missed in the 
numerical results of ^26. . 

For /3 = 0.49, the local bifurcation boundary of the 2nd family occurs at s = 0.894, and the eigenvalues A = 
±i [(3i{s) — /?2(s)] are isolated, such that vector solitons near the local bifurcation boundary are spectrally 
stable. This is indeed confirmed in Fig. 6, where the spectrum diagram is displayed for all values of s. It is 
seen that at the local bifurcation boundary s — 0.894, there are two pairs of isolated imaginary eigenvalues. 
The pair A = ±z [(3i{s) — /32(s)] = ±0.292z has the negative energy, while the pair A = ±0.310z has the 
positive energy. The second pair corresponds to the positive eigenvalue of £s in H4.9|l . As s moves leftward 
from the boundary point 0.894, these two imaginary eigenvalues move toward each other. At s = 0.891, 
they coalesce and create a quadruple of complex eigenvalues, in agreement with Appendix A. 3. However, 



3, f3ilL<P<h 
2, (3i<P<Pill, 



1, f3ilL<f3<l. 
0, Pi<f3- 



^ -"stab- 



(4.11) 



(4.12) 
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this instability persists only in a tiny interval 0.882 < s < 0.891, and it is very weak, with growth rates 
below 0.01. At s = 0.882, these complex eigenvalues coalesce and bifurcate back into two pairs of purely 
imaginary eigenvalues again. When s decreases further, one pair of these imaginary eigenvalues (denoted as 
±(Ti in Fig. 6) always remain imaginary, but the other pair of imaginary eigenvalues (±(72 in Fig. 6) move 
toward zero and become real for s < 0.80, in agreement with Appendix A.l. There is one more eigenvalue 
(0-3) in Fig. 6, which bifurcates from the edge of the continuous spectrum at s = 0.718, and always stays 
imaginary. The pattern of Fig. 6 differs from that of Fig. 5 in that complex instability is not set in at 
the local bifurcation boundary P — /32(s), and, once it is set in, it is confined in a narrow interval of s. 
Similar to the /? = 0.16 case above, this narrow interval of oscillatory instability was missed in but the 
exponential instability was captured there. 

Finally, we have mapped out the regions of exponential and oscillatory instabilities in the entire domain 
of existence for the 2nd family of vector solitons. The results are shown in Fig. 7. The almost-straight 
boundary lines show the local and nonlocal bifurcation boundaries [5]. The large domain of exponential 
instability away from the local bifurcation boundary corresponds to the one computed in Fig. 2 of j25j . 
The domain of oscillatory instability consists of two sub-domains. The larger sub-domain is at s < 0.857, 
where the eigenvalues A = ±i [/3i(s) — ^2(3)] at the local bifurcation boundary (3 = (32{s) are embedded. 
The smaller sub-domain is at s > 0.857, where the eigenvalues A = ±i[(3i{s) — /32(s)] are isolated but 
bifurcate to complex eigenvalues away from the local bifurcation boundary /3 = (32{s). Both sub-domains 
of oscillatory instabily were missed in . 

5 Summary 

We have applied the Closure Theorem for the negative index of the linearized Hamiltonian to multi-hump 
vector solitons in the general coupled NLS equations H2.1() . Unstable eigenvalues of the linearized problem 
(|2.4|l are approximated with the perturbation series expansions and found numerically with the shooting 
method. Not only the numerical results are in excellent agreement with the closure relation H2.9|l . but also 
the closure relation (|2.9(l shows that all unstable eigenvalues are recovered with the numerical shooting 
method. These analytical and numerical results establish that all multi-hump vector solitons in the non- 
integrable coupled cubic NLS equations are linearly unstable, while multi-hump vector solitons in the 
coupled saturable NLS equations can be linearly stable in certain regions of the parameter space. In the 
latter case, we have also discovered a new oscillatory instability which was missed before. This oscillatory 
instability significantly reduces the stability domains of vector solitons. 

We note that the Closure Theorem is applied differently in Sections 3 and 4. For coupled cubic NLS 
equations, we compute the closure relation H2.9|l from the right-hand side n{Ci) — piU) + n{Co) at the 
local bifurcation boundary and then match it with the number of unstable eigenvalues of the linearized 
problem t\'2A\\ . For coupled saturable NLS equations, we compute the closure relation 1)2. 9|l directly from 
the left-hand side N^eai + SA'comp + '^■^haag ^'^ the local bifurcation boundary and continue it in the entire 
existence domain. 

While the negative index theory is well understood in [JT] and well illustrated in ^1 E| and this paper, 
the following problem still remains a challenge: how do we understand the stability of vector solitons in 
sign-indefinite coupled-mode equations, such as the nonlinear Dirac equations? The Closure Theorem is 
obviously invalid for the Dirac equations as the continuous spectrum has both positive and negative energies. 
Thus a generalization of the Closure Theorem to such systems is highly desirable for future advances. 

Acknowledgements. The work of D.E.P. was supported in part by NSERC grant No. 5-36694. The work 
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A Appendix: Bifurcations of unstable eigenvalues 

Here we classify the four special cases, when one of the assumptions (i)-(iv) of the Closure Theorem fails. 
We derive sufficient conditions, when new unstable eigenvalues with Re(A) > bifurcate in the linearized 
problem 1)2.4(1 from eigenvalues with Re(A) — 0. For clarity of notations, we use an equivalent form of the 
spectral problem 1(2.4(1 : 

£iu = -/C^^u, ueX^(")(R), (A.l) 
where 7 = — A^ and is the constrained subspace of L^(R): 

= {u e i2 . (^^^^^^^ =0,n = 1, N} . (A.2) 

Here ei,...,ejv are unit vectors in and none of the components ^nix) is assumed to vanish identically 
on X e R. Operator Cq is always invertible in xi"-'(R), since eigenvectors {$„(a;)e„}^^]^ form a basis in 
the kernel of Cq. In the domain = {X E C : |A| > e} for any e > 0, there exists a relation between u(x) 
and w(x): 

w = A£o ^u, A e V„ (A.3) 

such that two eigenvectors (u, w)-^ and (u, — w)-^ of the linearized problem ((2.411 for A and (—A) correspond 
to a single eigenvector u(a;) of the problem IIA.l(l for 7 = — A^. 



A.l Zero eigenvalues of Ci and £0 

The kernel of Cq has a basis of N eigenvectors {<I>„(a;)e„}^^]^. Therefore, the existence domain of the 
vector solitons ((2.2(1 is confined by the boundaries, where <i>„(a;) = for some n = 1,...,A^. Let B be an 
open simply-connected domain in the parameter space /3jv), where the vector solitons ((2.2(1 exist. 

No bifurcations of zero eigenvalues of Lq occur in f3 & B. 

The kernel of Ci has always the eigenvector ^'(x). When z{Ci) — 1, this is the only eigenvector in the 
kernel of Ci. When z{Ci) > 1, additional linearly independent eigenvectors Uo(a;) exist in the kernel of Ci, 
such that the geometric multiplicity of A = in the linearization problem ((2.4(1 exceeds {N + 1). Under 
parameter continuation, the zero eigenvalue A = generally moves either to the real or purely imaginary 
axes of A. We study this bifurcation in the case, when z{£i) = 2, z{L() = 0, and (3 G B. 

Proposition A.l Let e be the bifurcation parameter and, at e — 0, there exists a non-zero eigenvector 
uo e X^'^\R), such that CiUq = and Uq is linearly independent of^'{x). Assume that Ci{e) and Co{£) 
are C^-functions at e — 0, such that Iq — (uo, £,^^(0)uo) ^ and Sli = (uo,£i(0)uo) 7^ 0. Then, there 
exists eo > such that the linearized problem {2.4}) has a real positive eigenvalue A in the domain: 

= {e : < |e| < eq, sign(e) = -sign(/o^Zi)} . 

Proof. We expand solutions of l(A.l(l in power series of e: 

u(x) -uo(a:) + eui(a:)+0(e2), j^eji+Oie^). (A.4) 
The function Ui (x) solves the non- homogeneous problem in Xc"'' (R) : 

A(0)ui+/:l(0)uo =7i'Co'(0)uo. (A.5) 
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Using the Fredholm Alternative Theorem, we find from (|A.5p that Sh = 71/0- If Sli ^ 0, la 0, and 
sign(e) = — sign(/o(5Zi), the eigenvalue 7 is negative in the first order of e, such that A = ±^—7 are real. ■ 

Corollary A.l Let Cq be positive definite, such that Iq > 0. A new negative eigenvalue /i(e) of Ci as e ^ 
results in a new negative eigenvalue 7(e) of the problem l^A.l]} . such that 

lim 1 = ^ ^ . (A.6) 

Bifurcation z{Ci) > 1 may occur only if iV > 1 in the system H2.1|l . Analysis of this bifurcation with the 
Lyapunov-Schmidt reduction method are reported in a similar content in |17l 1181 IHH] . 

A. 2 Zero eigenvalues of U 

When z{U) > 0, algebraic multiplicity of A — exceeds (2A^ + 2). Under a parameter continuation, the zero 
eigenvalue A = generally moves either to the real or purely imaginary axis of A. We study this bifurcation 
in the case, when z{U) = 1, z(£i) — 1, and f3 E B. 

Proposition A. 2 Let e be the bifurcation parameter and, at e = 0, there exists a non-zero eigenvector 
i> e M^, Uu = Q, such that the eigenvector uq G Ari"''(M) solves the problem: 

N 

CixiQ = - ^ j/„$„(a;)e„. (A.7) 

Assume that Ci{e), CQ{e), andU{e) are C^-functions at e ~ 0, such that Iq = (uq, £q ""^(Ojuo) 7^ and 
Su — {u,U'{Q)u) 7^ 0. Then, there exists eq > such that the linearized problem {2.4}) has a real positive 
eigenvalue A in the domain: 

= {e : < |e| < eo, sign(e) = -sign(Zo<^w)} • 

Proof. If there exists v E M^, such that U{0)i' — 0, then the eigenvector Uo{x) for the problem (|A.7p is 
given explicitly as 

, , d^(x) , . 

uo(x) = p. (A.8) 



dp.. 

Using Cid^ /dPn = —^n^n and Co^n^ti = for any e, we obtain the following derivative relations 



£o(0)*;(0)e,. +£i(0)4>„e,. = 0, (A.IO) 

where $'(0) stands for derivative of $(e) in e. We expand solutions of ljA.l(l in power series of e: 

u(x) =uo(a;) + eui(x)+0(e2), ^^^^^ + 0{e^). (A.ll) 

Since the eigenvector uo(x) solves the non-homogeneous problem l)A.7|l . the relation between u(a;) and w(a;) 
is modified as follows: 



1 ^ 

w = A£o^u+ -^z/„$„(x)e„. (A.12) 



n=l 
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As a result, the function Ui(a;) solves the non-homogeneous problem: 



N 



£i{0)ui + £'i(0)uo = 7i£o i(0)uo - ^ i.„<(0)e„, 



r! = l 



subject to the constraints: 

($„(0)e„,ui) + «(0)e„,uo) - 0. 
Using the Fredholm Alternative Theorem, we find from (|A.13|1 and (|A.14|I that 



N 



(uo,/:'i(0)uo) =7i(uo,/:o"'(0)uo) -2^z.„«(0)e„,Uo> 



11=1 



As a result, 



N N 



7l(Uo|£o (O)Uo) = 2^ 2^ i^nl^m^(*m(e)em, 
n— 1 m— 1 



(j^,Z^'(0)i/) 



(A.13) 
(A.14) 

(A.15) 
(A.16) 



e=0 



such that Su — 2jiIq. When Su 7^ 0, Zo 7^ 0, and sign(e) — — sign(Zo<5u), the eigenvalue 7 is negative in the 
first order of e, such that A = ±-^—7 are real. ■ 

Corollary A. 2 Let Co be positive definite, such that lo > 0- A new negative eigenvalue ii{e) of hi as e ^ Q 
results in a new negative eigenvalue 7(e) of the problem \A.l\l . 

Bifurcation z{U) > was analyzed in |2HII2] with power series expansions of the problem (|2.4I) near A = 
and in P rro' with Taylor series expansions of the Evans function. 



A. 3 Multiple non-zero eigenvalues of zero energy 

When the problem 12.4|l has a multiple eigenvalue A = Ao G iM of zero energy, the corresponding eigenvector 
(uo,wo)"^ satisfies the conditions of the Fredholm Alternative Theorem: (uo,£iUo) = 0, (wo,>CoWo) = 0, 
such that ft,[uo,Wo] = and lo — (uq, >C|^"'^Uo) = 0. Under parameter continuations, multiple eigenvalues 
are generally destroyed and new complex eigenvalues A may arise in the problem H2.4() . We study this 
bifurcation in the case, when a multiple eigenvalue A = Aq has algebraic multiplicity two and geometric 
multiplicity one, while /3 e S U dB. 

Proposition A. 3 Let e be the bifurcation parameter and, at e = 0, there exist non-zero eigenvectors 
uo, Ug e Xc'^\r) for 70 e M and 70 < /?^i„, such that 

£iUo = 7o£(7^uo, (A. 17) 

Liu'o = jn£o^u'„+Co^un, (A.18) 

and lo = (uo,>C^^Uo) — 0. Assume that Ci{e) and Co{e) are C^-functions at e = 0, such that — 
(uo, £(7^(0)uq) ^ and Sh ~ (uq, {C[{0) — jo^o^'i^)) "o) 7^ 0. Then, there exists eo > such that the 
linearized problem ^2.4^ has two complex eigenvalues A with Re(A) > in the domain: 

= {e : < |e| < eq, sign(e) = ~sign{l'o6h)} . 

Proof. We expand solutions of ljA.l|) in power series of e^/^: 

u{x) = Uo{x)+e^/^-fiu'o{x)+€U2{x) + 0{e^/^), (A.19) 
7 = 70 + e^^^7i+ £72 + 0(6^/2). (A.20) 
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The function U2 (a;) solves the non- homogeneous problem in Xc' (R) : 

/:i(o)u2 + /:'i(o)uo -7o/:^'(o)u2 + 7o/:o-i'(o)uo + 7?/:o"'(o)u;, + 72/:o'(o)uo. (A.21) 

Using the Fredholm Alternative Theorem, we find from ljA.2ip that 6h = 71^0. Since (7 — 70)^ ~ e7i + 
0(e'^/^), the eigenvalues 7 bifurcate into the complex plane if /q ^ 0, Sh ^ 0, and sign(e) = — sign(/Q(5ft,). ■ 

Corollary A. 3 Let Iq ^ Q and Sh ^ 0. There exists eo > such that the problem I^A.l]) has two real 
eigenvalues 7 in the domain, 

V, = {e : < \e\ < eo, sign(e) = sign(Z^(5/i)} 
with oppositely signed quadratic forms 

(u,£o-iu) =2ei/VZo + 0(e). 

When < 70 < Pmin, multiple eigenvalue A = Ao is purely imaginary, and the bifurcation of Proposition^^! 
is the instability bifurcation. When 70 < 0, the eigenvalue A = Aq is purely real and is thus already unstable. 
Characteristic features of bifurcation of multiple eigenvalues were analyzed in |12) . This bifurcation is 
generic when purely imaginary eigenvalues of positive and negative energies /i[u,w] coalesce, according to 
Corollary I A. 31 1321 . General results on the collisions of purely imaginary eigenvalues of different energies 
and the concept of so-called Krein signatures can be found in |22| ■ 

A. 4 Embedded eigenvalues 

When the problem 12.411 has an embedded eigenvalue A = Aq £ «IR with |Im(Ao)| > Pmin, it is generically 
unstable under parameter continuation EUJ |^ . When it has a positive energy, it disappears from the 
continuous spectrum, while when it has a negative energy, it bifurcates as complex unstable eigenvalues 
with Re(A) > jH]- We study this bifurcation in the case, when an embedded eigenvalue has the geometric 
and algebraic multiplicities one, while /3 E BU dB. 



Proposition A. 4 Let e he the bifurcation parameter and, at e — 0, there exist a non-zero eigenvector 
uo G xi"^(M) for 70 G M, 70 > /3^in, such that 

Ciuo = 7o/:o 'uo. (A.22) 

Assume that Ci{e) and Cq{£) are -functions at e ~ 0. such that Iq ~ (ug, £q ^ (0)uo) < and T ^ in 
I^A . 29\} below. Then, there exists eq > such that the linearized problem \2.4\l has two complex eigenvalues 
A with Re(A) > 0. 

Proof. For embedded eigenvalues, we use the linearized problem in the original form p. 4(1 . Consider 
perturbation series expansions near the embedded eigenvalue A = Ao, with Im(Ao) > /3min: 

nix) = Uo(.x) +eui(a;) +e2u2(a;) + 0(e3), (A.23) 

w(a;) = wo(a;) + ewi(x) + eV2(a:) + 0(e^), (A.24) 

and 

A = Ao + eAi+e2A2 + 0(e3). (A.25) 
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Corrections of the perturbation series (|A.23p and (|A.24|I satisfy linear non-homogeneous equations following 
from the linearized problem H2.4|l : 

£i(0)ui + Aowi = -£'^(0)uo - AiWo, 

/:o(0)wi-AoUi = -4(0)wo + AiUo (A.26) 

and 

£i(0)u2 + A0W2 = -£i(0)ui - i£"(0)uo - Aiwi - A2W0, 

£o(o)w2-AoU2 = -/:^,(o)wi-i/:;;(o)wo + Aiui + A2Uo. (A.27) 

Bounded solutions of the problem (|A.26|I exist only if the right-hand-side is orthogonal to the eigenvector 
(uo,wo). The solvability condition results in the equation: 

Ai ((wo,uo> - (uo,wo>) = (uo,£i(0)uo> + (wq, £o(0)wo), 

such that Re(Ai) = 0. Since the eigenvalue A = Ao belongs to the continuous spectrum of the problem (|2.4|l . 
the correction terms (ui, Wi) have non-vanishing tails in the limit \x\ — > oo. Assuming that f3i < P2 ^ ■■■ < 
Pn, we add Somnierfeld radiation conditions to uniquely determine the correction terms (ui, wi): 



where gj are some constants, kj = (Im( Aq ) — Pj)/ dj , and Kx^ is the number of branches with kj e 
It follows from (IA.26I) that 



lm((ui,£i(0)uo) + (wi,£o(0)wo) + Ai(ui,Wo) - Ai(viri,uo)) 
((ui,£iui) + (wi,£owi) - (£iui,ui) - (£owi,wi)) 



if AO 

= -2^d,fc,(|g+|2 + |5T|2)^r<0. (A.29) 

Again, bounded solutions of the problem ljA.27|l exist only if 

A2 ((wo, uo) - (uo, v^fo)) + Ai {(-wq, Ui) - (uq, wi)) 
= (uo,£'i(0)ui) + (wo,£[,(0)wi> + i(uo,£'/(0)uo) + i(wo, £;,'(0)wo), (A.30) 

such that 

^^^^^^ ^ Im((ui,A(0)uo) + (wi,i:;,(0)wo)) j^ d,k,{\g+\^ + \g-\^) ^ 
^ 2Im(Ao)(uo,£o ^uo) ^ Im(Ao)(uo, ^q^^Uo) 

When r 7^ and Iq = {uq, Cq^{0)uo) < 0, the embedded eigenvalue A = Aq becomes a complex unstable 
eigenvalue A with Re(A) > 0. ■ 

Corollary A. 4 The linearized problem {2.4^ does not have complex or embedded eigenvalues X if lo — 
(uo,/:o ^(0)uo) >0 and r^O. 

Proof. The formal computation in (|A.3ip predicts that Re(A2) < 0, when F ^ and > 0. However, the 
correction terms (ui,Wi)'^ in (jA.28l) grow exponentially in x, since fcj(A) = ^(Im(A) -iRe(A) -/3j)/dj 
implies that Im(fcj) > 0. The embedded eigenvalue A = Aq becomes a resonant pole with Rc(A) < 0. ■ 
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Characteristic features of bifurcations of embedded eigenvalues were analyzed in |H1 El ■ This bifurcation 
is generic for multi-hump vector solitons at the boundaries of the existence domain (3 S dB |35| . 



Summarizing, there exist four bifurcations, which may lead to unstable eigenvalues in the spectral problem 
(|2.4|l : (i) z{Ci) > 1, (u) ziU) > 0, (in) multiple eigenvalue Ao £ iM, |Im(Ao)| < /3min, and (iv) embedded 
eigenvalue Aq S «M, |Im(Ao)| > Pmin- Let nx{h) be the negative index of the linearized Hamiltonian in the 
constrained space Xc"''(M), i.e. the number of negative eigenvalues of Ci and £o in xi^^{R). It is known 
from [ni (see also |13EH1) that 

nx{h)^n{Li)-p{U)+n{Co). (A.32) 

The closure relation ()2.9|l gives then: 

nx (h) = A^rcai + ^N,-^^^ + 2A^eomp. (A.33) 

It is clear from (|A.33(I that bifurcations (i) and (ii) change the negative index nx{h) due to a change in 
Nrcai by one, while bifurcations (iii) and (iv) do not change the negative index nx{h) due to an exchange 
in ^imag and TVcomp. 
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Figure 1: Vector solitons of the Ist-family in the cubic model (|3.1|) a± (3 — 0.36 and three x values marked 
as "*" in the upper left figure. Here uo = y^, and {u,v) = (<i>i, $2)- 




Figure 2: Eigenvalue spectrum of the 1st family of vector solitons in the cubic model (|3.1(l at /? = 0.36. 
The solid (dashed) lines are the real (imaginary) parts of the eigenvalues. 
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Figure 3: Vector solitons of the 2nd family in the cubic model (|3.1() at /3 = 0.36 and three x values. 
Notations ui, u and v are the same as in Fig. 1. 
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Figure 4: Eigenvalue spectrum of the 2nd family of vector solitons in the cubic model H3.1|l at /? = 0.36. 
The solid (dashed) lines are the real (imaginary) parts of the eigenvalues. 
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Figure 5: Eigenvalue spectrum of the 2nd family of vector solitons in the saturable model (|4.1|l at /3 = 0.16. 
The solid (dashed) lines are the real (imaginary) parts of the eigenvalues. 




Figure 6: Eigenvalue spectrum of the 2nd family of vector solitons in the saturable model (|4.1|l at /3 ~ 0.49. 
The solid (dashed) lines are the real (imaginary) parts of the eigenvalues. The right figure is a zoom-in of 
the left figure. 
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Figure 7: Regions (shaded) of exponential and oscillatory instabilities of the 2nd family of vector solitons 
in the saturable model ()4.1|l . 



24 



